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SUMMARY 


The  problem  of  testing  whether  one  regression  function 
is  larger  than  another  on  a  specified  compact  set  R  is  consi¬ 
dered.  The  regression  functions  must  be  linear  functions  of 
the  parameters  but  need  not  be  linear  functions  of  the  inde¬ 
pendent  variables.  The  proposed  test  statistic  is  compared  to 
a  standard  t  percentile.  The  test  has  an  exactly  specified  size  in 
typical  situations.  Properties  of  the  power  function  of  the 
test  are  investigated.  The  related  question  of  comparing  a 

regression  function  to  a  specified  function  is  also  considered.  j 


By- 
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1.  Introduction 


! 
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In  this  paper  we  consider  testing  whether  the  regression 
function  from  one  population  is  everywhere  above  the  regression 
function  from  another  population.  A  medical  researcher,  for 
example,  might  be  interested  in  testing  whether  the  mean  of 
a  response  variable  in  a  diseased  population  is  larger  than  the 
mean  of  a  response  variable  in  a  healthy  population  for  all  pos¬ 
sible  values  of  an  independent  variable.  Examples  such  as  these 
have  been  discussed  by  Tsutakawa  and  Hewett  (1978),  Hewett  and 
Lababidi  (1980)  and  Spurrier,  Hewett  and  Lababidi  (1980)  who  have 
considered  this  testing  problem. 

The  model  considered  herein  generalizes  the  models  of 
Tsutakawa  and  Hewett  (1978)  and  Hewett  and  Lababidi  (1980)  in 
two  ways.  First  the  regression  functions  may  be  functions  other 
than  linear  functions  of  the  independent  variables.  For  example, 
the  regression  functions  may  be  quadratic  or  higher  degree  poly¬ 
nomials.  Second,  the  two  regression  functions  need  not  be  of 
the  same  functional  form.  For  example,  one  may  be  assumed  to 
be  a  linear  function  and  the  other  a  quadratic  fun.  The 

test  proposed  herein  reduces  to  the  tests  of  Tsutakawa  and 
Hewett  (1978)  and  Hewett  and  Lababidi  (1980)  for  the  special 
models  they  consider.  The  proposed  test  requires  no  new  tables 
for  its  implementation.  Only  a  standard  t  table  is  needed. 
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The  model  and  test  are  presented  in  Section  2.1.  A 
numerical  example  using  the  organ  weight  data  of  normal  and 
diabetic  mice  is  presented  in  Section  2.2.  Properties  of  the 
test,  including  its  power,  are  discussed  in  Section  3.  The 
related  problem  of  testing  whether  a  single  regression  func¬ 
tion  everywhere  exceeds  a  specified  function  is  discussed  in 
Section  4. 

2.  Model,  Test  and  Application 
2.1.  Model  and  Test 


Let  { CXj . ,  ) ,  j  *>  1,  ....  nj)  and  {(X^,  Y^);  j  *  1,  ...,  n2> 

denote  two  independent  sets  of  observations  where 
X. .  =  (X. ,  ....  X. .. ) .  The  X. .  may  be  observed  random  vectors 
or  design  variables  fixed  by  the  experimenter.  The  entire  analysis 
is  conditioned  on  the  observed  values  of  X^^ .  Let  R  denote  a 
closed  and  bounded  subset  of  k-dimensional  Euclidean  space.  R 
is  the  set  of  possible  values  of  the  independent  variables  X^. 

We  assume  that  given  the  X„  the  Y^  are  independent  normal  ran¬ 
dom  variables  with 


E(Y.  .  |X.  ..  =  x.  ,  ....  X. .  »  x. ..  ) 

v  13'  13I  13I  ’  13k  13k' 


l.*im  fim(xijr  *  *  *  *  xijkJ  “  fi(xij)6i 


m*l 


and 


var(,ijixiji  ■  xm . V  *  *ijk>  ■ 0  • 

2 

0.  )',  i  •  1,  2,  and  o  are  unknown  parameters. 

lpi 


1*. 


0j*(8£2  •  •  •  • » 


•  •  •  » 
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The  f^x)  =  (f.  (x),  ....  f.  (x)),  i  =  1,  2,  are  known  vectors 

"i 

of  functions  which  define  the  functional  form  of  the  regression 

functions.  For  example  the  f^j  might  be  polynomials  such  as 
2 

1,  Xj,  x2,  or  XjXj.  By  allowing  px  *  p2  and  flm  *  f2m(x), 
this  model  allows  the  two  regression  functions  to  have  different 
functional  forms.  The  first  might  be  a  linear  function  and  the 
second  a  quadratic  function.  But  usually  Pj  *  P2  and  fj(x)=  f2(x) 
will  be  chosen  so  the  regression  functions  have  the  same  func¬ 
tional  form. 

We  wish  to  compare  the  regression  functions  f1(x)B1  and 
f 2(x)B2-  In  particular  we  are  interested  in  whether  fjCxiBj 
is  always  greater  than  f2(x)B2.  The  test  we  will  propose  is 
a  size  a  test  of 

Hq:  f^xjBj  s  ^or  at  least  one  x  €  R 

vs. 


V  fjW&j  >  f2(x)e2  for  every  x  c  R. 

Let  bj  and  b2  denote  the  least  squares  estimates  of  ^  and 


B2  and  let 


2  n* 


*  •  I  l  (Xij  -  fi(*i i)V  /v 
2 

denote  the  pooled  estimate  of  a  where  v  ■  nj  "  Pj  *  n2  "  P2* 

The  estimate  b^  has  a  multivariate  normal  distribution  with  mean 

2  -1 

Bj  and  covariance  matrix  o  where  the  (m,  n)  element  of  D.  is 
ni 

WW'  ut  'W  -  fi«»;  *  qooo;1  qoo. 


I 


Then  the  variance  of  fjMbj  -  £2 (x)b2  is  c  e(x). 

The  test  we  will  propose  may  be  motivated  in  this  way. 
Consider  comparing  the  two  regression  functions  only  at  a  sin¬ 
gle  point  x  e  R.  Let  Tx  *  (fjMbj  -  f2Cx)b2)/s»/,eCx)  •  The 
test  which  rejects  HQx:  f^xJBj  s  f2(x)02  in  favor  of 

HAx:  *j(*)Bj  >  When  Tx  >  ti  a  (v)  is  a  s^ze  °  test 

where  t.  (v)  is  the  (1  -  o)th  percentile  of  a  t  distribution 
l-a 

with  v  degrees  of  freedom.  The  alternative  of  interest  H  is 

the  intersection  of  all  the  H.  for  all  x  c  R.  It  seems  rea- 

Ax 

sonable  to  decide  in  favor  of  only  if  the  test  bnsed  on 
Tx  decides  in  favor  of  for  every  x  e  R.  This  leads  us  to 

propose  the  following  test. 

Define  the  test  statistic  T  by 

(2.1)  T  =  min  T  . 

xeR  * 

Reject  Hq  in  favor  of  HA  if  and  only  if  T  >  ^  0(v)*  This  test 
is  always  a  level  a  test  in  that  the  probability  of  a  type  one 
error  is  always  less  than  or  equal  to  a.  This  test  has  size 
exactly  equal  to  2.  if  the  f^  are  continuous  functions  of  x 
and  there  are  values  of  Bj  and  B2  such  that  fj(x)8j  »  f2(x)B2 
for  one  value  of  x  e  R  and  fjCxjBj  >  f2(x)B2  for  all  other 
x  €  R.  These  facts  are  proven  in  the  Appendix.  Section  3  con- 
tails  examples  of  when  the  condition  is  satisfied  and  the  size 
is  exactly  a. 
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The  test  statistic  T  will  usually  have  to  be  evaluated 

by  numerical  methods.  This  is  discussed  in  Section  3.  In 

some  cases  it  can  be  evaluated  explicitly.  In  particular  it 

is  shown  in  the  Appendix  that  this  test  is  equivalent  to  the 

tests  proposed  by  Tsutakawa  and  Hewett  (1978)  and  Hewett  and 

Lababidi  (1980)  for  the  special  models  they  consider  if  a  s  .5. 

2.2.  A  Numerical  Application 

Ne  consider  the  data  relating  body  and  kidney  weight  for 

diabetic  and  healthy  mice  which  was  presented  by  Bishop  (1973) . 

These  data  were  analyzed  by  Tsutakawa  and  Hewett  (1978) .  The 

data  are  plotted  in  Figure  1.  The  independent  variable  x  is 

the  body  weight  and  the  dependent  variable  y  is  the  kidney 

weight.  We  consider  modelling  the  mean  of  the  kidney  weight  as 

a  quadratic  function  of  the  body  weight.  We  consider  testing 

2  2 

H0:  B11  *  612X  *  B13x  s  821  *  B22X  *  B23X  for  some  x  ‘  R 

2  2 

versus  HA:  8n  ♦  B12x  ♦  BJ3x  >  &21  ♦  B22x  ♦  e23x  for  all  x  c  R 

where  the  first  population  is  the  diabetic  population  and 

R  «  (x:  26  <  x  s  52}.  The  least  squares  line  for  the  9  dia- 

2 

betic  mice  is  given  by  1821.14  -  49.92x  ♦  .77x  .  The  least 
squares  line  for  the  25  healthy  mice  is  given  by 
-434.68  ♦  62.70x  -  .85x2.  s2  =>  (66052  ♦  196505)/(9  -  3  ♦  25  -  3) 
«  9377.0.  The  matrices  Dj1  and  D*1  are 
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and 


408.0894 

-19.1506 

.2214 

-19.1506 

.9031 

-.0105 

.2214 

-.0105 

.0001_ 

"81.1483 

-4.6954 

.0666 

-4.6954 

.2735 

-.0039 

.0666 

-.0039 

.0001 

Thus  Tx  *  a(x)/s/e(x)  where  a(x)  =  2255.82  -  112. 62x  ♦  1.62x2 
and  e(x)  =  489.2377  -  47.6920x  ♦  l.7526x2  -  .0288x3  ♦  .0002x4. 
Numerical  minimization  of  for  x  between  26  and  52  yields  the 
test  statistic  T  *  .4797  which  corresponds  to  x  »  40.10.  This 
T  is  at  approximately  the  68th  percentile  of  a  t  distribution 
with  28  degrees  of  freedom. 

3.  Properties  of  the  Test 
3.1  Power  Function  Properties 

Tho  test  we  propose  has  size  exactly  o  if  the  f^  are  con¬ 
tinuous  functions  and  there  are  values  of  and  62  such  that 

fl^6l  *  *2^®2  f°T  °ne  value  x  e  R  80(1  fjfxJgj  >  f 2 (x) 02 
for  all  other  x  c  R.  This  condition  is  satisfied  if  the  f^(x) 

include  the  constant  1,  the  linear  functions  x.,  i  *  1,  ...,  k, 

and  the  quadratic  functions  x^x^ ,  i  a  i,  . k,  j  a  1,  ...,  i. 


Then  gj  and  g2  can  be  chosen  so  that 

*  V<x  "  x0^’  the  s<luare  of  the  <iis' 

tance  from  x  to  x^,  where  x^  is  a  fixed  element  of  R.  For  this 

choice,  fj(x)gj  -  f2(x)g2  is  zero  for  x  =  x^  and  is  positive  for 

all  other  x.  Another  situation  in  which  the  condition  is  satis- 

k 

fied  is  if  f.(x)g.  3  Bi(k+1)  ♦  and 

R  »  {x:  x.  <  x.  £  x  *  j  =  1,  ....  k},  the  model  considered 

J  *  J  J 

by  Tsutakawa  and  Hewett  (1978)  and  Hewett  and  Lababidi  (1980). 


The  choice  of  B2  3  0,  g^ 


=  8lk  *  1  and  6l(k+l)  = 
c.  -  x.J  which  is  zero  for 


yields  fJ(x)g1  -  f2(x)g2  *  £  (x^  -  x^)  which  is  zero  for 
x  3  Cxj*»  •  xk*)  and  positive  for  all  other  x  e  R.  A  numeri¬ 

cal  corroboration  of  the  fact  that  the  size  is  exactly  a  is 
found  in  the  following  simulation  study. 

A  simulation  study  was  conducted  to  investigate  the  power 

function  of  the  test.  In  this  study  the  regression  functions 

2  2 
were  fj(x)g^  *  g^  ♦  g^2x  +  3^x  ,  is  i,  2.  The  variance  o 

was  set  equal  to  one.  The  sample  sizes  n^  and  n2  were  both  10 
with  3  observations  at  each  of  x  =  1  and  x  3  -1  and  4  observa¬ 
tions  at  x  »  0.  R  s  (x:  -1  s  x  <  l).  The  size  of  the  test 
was  fixed  at  as  .05  by  using  t  95(14)  *  1.761.  The  International 
Mathematics  and  Statistics  Library  programs  GGNSM  and  GGCHS 

were  used  to  generate  the  random  vector  b^  -  b2  and  the  random 
2 

variable  S  .  A  total  of  3000  repetitions  were  used  to  obtain 
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each  of  the  estimates  in  Tables  1  and  2. 

The  maximum  probability  of  a  Type  I  error  takes  place 
when  fjtxJgj  *  f2(x)62  for  one  x  and  f^x^  -  f2(x)82  becomes 
large  for  all  other  x.  This  can  be  observed  in  Table  1  where 
the  probability  of  a  Type  I  error  is  given  for  various  values 
of  and  *n  t*ie  null  hypothesis.  As  one  proceeds  down 
Columns  II,  III,  IV  or  V  of  Table  1,  fj(x)8j  «  f2(x)&2  ^or  one 
value  of  x(x  =  -1  for  Columns  III,  IV  and  V  and  x  =  0  for 
Columns  II)  and  f^(x)g^  -  f2(x)82  is  becoming  large  for  all 
other  values  of  x.  The  probability  of  a  Type  I  error  increases 
to  a  a  .OS  as  one  proceeds  down  any  column.  The  estimates 
slightly  exceed  .05  in  a  few  cases  due  to  sampling  error. 

The  power  function  of  the  proposed  test  exhibits  the 
following  monotonicity  property.  If  (8^,  B2)  and  (Bp  8*)  are 
two  parameter  vectors  which  satisfy 

(3.1)  f1(x)8j  -  f2(x)8j  *  fjCxJfij  -  f,(x)B2 

for  every  x  with  strict  inequality  for  some  x,  then  the  power 
at  (Bp  Bp  is  greater  than  the  power  at  (Bji  &2) •  This  pro¬ 
perty  is  apparent  as  one  proceeds  down  any  column  in  Tables 
1  or  2,  across  the  first  three  columns  in  any  row  of  Table  2, 
or  across  the  last  three  columns  in  any  row  of  Table  1.  It  is 
also  apparent  if  one  compares  any  two  corresponding  entries  in 
Table  2a  and  2b  since  the  regression  functions  are  farther  a- 
part  in  2b  than  in  2a. 
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The  power  of  the  test  is  near  one  only  if  fjCx)^  -  f2(x)82 
is  large  for  all  x.  This  is  the  case  for  the  lower  entries  in 
Table  2.  In  Table  2,  the  minimum  distance  between  the  regression 
functions  is  c.  The  power  nears  one  only  as  the  minimum  dis¬ 
tance  c  becomes  large. 

The  test  we  propose  is  biased  in  that  the  probability  of 
rejecting  H.  is  less  than  a  for  some  (8,,  8-)  in  H  .  The  fea- 
ture  was  noted  by  Tsutakawa  and  Hewett  (1978)  for  the  special 
model  they  considered  and  it  continues  to  exist  for  the  more 
general  models  we  consider.  This  biasedness  can  be  observed 
in  the  entries  for  c  =  .5  which  are  less  than  .95  in  Table  2. 

But  as  noted  by  Tsutakawa  and  Hewett  (1978)  for  their  special 
case,  the  test  we  propose  is  consistent  in  that,  for  any  fixed 
point  (8  ,  e2)  in  Na»  the  power  can  be  made  arbitrarily  near 
one  by  choosing  the  sample  sizes  sufficiently  large.  Although 
we  do  not  feel  this  bias  is  serious,  it  should  be  noted  that 
the  power  of  the  test  we  propose  may  be  small  if  fjWBj  ex¬ 
ceeds  f2(x)82  by  only  a  small  amount  over  most  of  R. 

The  power  function  properties  we  have  described  in  this 
section  are  true  in  general,  not  just  for  the  case  of  quadratic 
regression  we  considered  in  the  simulation  experiment.  The 
proofs  of  these  facts  can  be  accomplished  using  the  methods  em¬ 
ployed  in  proofs  in  the  Appendix. 


-  10  - 


3.2.  Computational  Shortcuts 

In  (2.1)  the  test  statistic  T  was  defined  as  the  minimum 
of  T^  over  the  set  R.  Typically  the  computation  of  the  test 
statistic  will  be  accomplished  by  a  numerical  minimization  of 
T  .  But  to  perform  the  test  the  actual  value  of  T  need  not 
be  computed.  One  only  needs  to  know  whether  T  >  t^  ^(v) 
or  T  i  tj  a(v).  In  this  section  we  describe  two  shortcuts 
which  allow  the  determination  of  whether  T  >  tlQ(v)  or 
T  s  tj  ^(v)  without  the  actual  computation  of  T. 

3.2.1.  Shortcut  for  determining  if  HQ  is  accepted. 

Let  X*  denote  an  arbitrary  finite  subset  of  R.  For  example, 
if  R  =  {x:  x ^  s  x^  £  x*,  i  =  1,  ...,  k),  X*  might  be  the  set 
of  2k  extreme  points  (x, ,  ....  x.)  where  x.  "  x.  or  x.  =  x?  . 

1  K  i  i  *  i  i 

Let  T'  *  min  T  .  Since  T  <  T',  if  T'  £  t.  (v)  accept  H„. 

xeX*  x  I"0  0 

Furthermore  if  T'  »  t'  the  significance  probability  associated 

with  T  is  at  least  P(TQ  >  t')  where  has  a  central  t  distri¬ 
bution  with  v  degree*  of  freedom. 

For  the  body  kidney  weight  data  of  Section  2.2, 

R  *  (x:  26  £  x  £  52).  T  for  x  *  26  is  .9962  and  T  for  x  =  52 

x  x 

is  .6347.  Either  one  of  these  points  is  less  than  T  g5(28)  =  1-7011 
so  the  test  accepts  HQ  at  level  .05.  Furthermore  the  significance 
probability  associated  with  T  is  at  least  P(T^  2  .6347)  *  .27. 

The  exact  significance  probability  computed  in  Section  2.2 


was  .32. 
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In  Tables  1,  and  2  the  second  (middle)  number  for  each 
entry  is  the  proportion  of  the  acceptances  which  were  deter¬ 
mined  by  the  shortcut  method.  These  values  indicate  that  the 


usefulness  of  this  shortcut  depends  on  the  actual  regression 
functions.  But,  in  many  cases,  a  large  proportion  of  the 
acceptances  were  determined  by  this  shortcut.  In  18  out  of 
the  62  cases  in  which  there  were  some  acceptances,  all  of 
the  acceptances  were  determined  by  the  shortcut. 

3.2.2.  Shortcut  for  determining  if  HQ  is  rejected. 

For  this  shortcut  to  be  valid, a  must  be  no  more  than  .5. 
Since  o  usually  satisfies  as  .1,  this  restriction  is  not  prac 
tically  important.  Let  m  denote  the  number  of  distinct  non- 
constant  functions  in  {^(x):  i  »  1,  2;  j  *  1,  ...,  p^) 

Let  Z*  denote  the  set  of  2m  points 


(zj,  z*)  »  C(zu,  ....  (*21>  •••»  *2p  formed  by 

replacing  f. ,(x)  by  either  max  f . . (x)  or  min  f . . (x)  in 
13  xeR  13  xeR  13 

(fj(x)»  f2(x)).  Note  that  if  fJr(x)  a  f2s(*)  then  zlr  s  Z2S» 

i.e.,  the  maximum  or  minimum  is  used  on  both  zj  and  z*.  Let 

T*  =*  min  T  *  where 
z*eZ*  z 


Tz*  “  (zlbl  '  Z2b2)/S  /*lDi1,r  *  Vl'Y  '  If  T*  >  Va(v) 

then  T  >  t.  (v)  and  HA  can  be  rejected.  Furthermore,  if 
i-a  0 

T*  «  t*,  the  significance  probability  associated  with  T  is  less 


-1  *- 


than  or  equal  to  P(Tq  >  t*)  where  TQ  has  a  central  T  distribu¬ 
tion  with  v  degrees  of  freedom.  If  T*  *  T where  z*  »  (f^Cx),  f2(x)) 
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for  some  x  e  R  then,  in  fact,  the  significance  probability 
equals  P(TQ  >  t*). 

2 

For  example,  suppose  fj(x)  *  (1,  x,  x  ),  and  f^x)  =  (1,  x), 

and  R  =  {x:  -1  s  x  <  2}.  There  are  two  distinct  nonconstant 

2  2 

functions,  x  and  x  ,  and  the  2=4  values  of  (z*,  zp  in  Z*  are 
(Cl.  -1,  0),  (1,  -1)),  ((1,  -1,  4),  (1,  -1)),  ((1,  2,  0),  (1,  2)) 
and  ((1,  2,  4),  (1,  2)).  Points  like  ((1,  -1,  0),  (1,  2)) 
where  x  has  been  replaced  by  its  minimum  in  z*  and  its  maximum 
in  t*  are  not  in  Z*. 

The  validity  of  this  shortcut  is  based  on  two  facts,  the 
fact  about  functions  which  are  the  ratios  of  linear  functions 
and  square  roots  of  positive  quadratic  functions  mentioned  in 
the  proof  of  Theorem  3  and  the  fact  that  A  »  {(fjfx),  f^x)):  x  e  R} 
is  a  subset  of 


B  a  {(2.,  z  ):  min  f..(x)  *  2..  *  max  f. .(x)  and  z  *  2  if 
L  *  xcR  13  13  xeR  3 

flr(x)  =  f2s(x)} 

so  a  minimum  over  A  is  not  less  than  a  minimum  over  B. 


This  shortcut  was  used  in  the  simulation  study  of  Section 

2 

3.1.  In  this  case  R  =  (x:  -1  s  x  s  1).  fj(x)  =  f2(x)  =  (1,  x,  x  ) 
so  there  are  m  »  2  distinct  nonconstant  functions.  The  2m  =  4 


points  in  Z*  are  ((1,  -1,  0),  (1,  -1,  0)),  ((1,  -1,  1),  (1,  -1,  1)), 
((1,  1,  0),  (1,  1,  0))  and  ((1,  l,  1),  (l,  i,  1)). 

In  Tables  1  and  2,  the  third  (bottom)  number  for  each  entry 
is  the  proportion  of  the  rejections  which  were  detected  by  this 
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shortcut.  The  usefulness  of  this  shortcut  is  seen  to  depend 
on  the  actual  value  of  the  regression  function  but  in  many  cases 
the  proportion  is  fairly  high.  In  17  of  71  cases  all  of  the  re¬ 
jections  were  detected  by  this  method,  avoiding  numerical  mini¬ 
mization  of  T  . 

x 

3.2.3.  Numerical  minimization  of  T  . 

x 

If  neither  of  the  shortcuts  presented  in  Sections  3.2.1  or 
3.2.2  determine  whether  is  to  be  accepted  or  rejected  or  if 
the  exact  value  of  T  is  desired  to  compute  an  exact  significance 
probability,  then  the  function  must  be  minimized  by  numerical 
methods  to  determine  the  value  of  the  test  statistic  T.  The  pro¬ 
blem  of  minimizing  a  function  such  as  which  is  the  ratio  of  two 
functions  of  x  has  been  studied  extensively  in  the  mathematical 
programming  literature  by  Charnes  and  Cooper  (1962) ,  Swarup  (1965) 
Sharma  (1967)  and  Craven  and  Mond  (1973,  1975a,  and  1975b).  These 
authors  have  found  that  this  non-linear  programming  problem  is 
equivalent  to  other  nonlinear  programming  problems  which  do  not 
involve  fractions.  These  results  could  simplify  the  numerical 

minimization  of  T  . 

x 

4.  Comparison  of  a  Regression  Function  with  a  Specified  Function. 

The  hypothesis  that  the  regression  function  from  a  popula¬ 
tion  is  everywhere  above  a  specified  function  can  be  tested  with 
a  test  similar  to  the  one  described  in  Section  2.  For  this  pro¬ 
blem  we  only  have  the  sample  {(X^,  Y^);  j  *  1.  . ...  from 


t. 
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the  first  population.  We  wish  to  test 

Hq:  ffxjflj  £  g(x)  for  at  least  one  x  e  R 
vs 

Ha:  fCxJfl^  g(x)  for  all  x  e  R 
where  g(x)  is  a  specified  function.  The  function  g(x)  may  be 
a  constant  c  in  which  case  we  are  testing  whether  the  mean  of 
the  response  is  greater  than  c  for  all  possible  values  of  the 
independent  variable  x.  Let 

nl 

’l  ■  -  WV2'”! 

where  Vj  ■  nj  -  pj4  Let  e1(x)  *  f^xJD^f'fx)  and 

TJx  *  (fjfxJbj  -  gCxJJ/Sj  /e^ (x) .  Define  the  test  statistic  T^ 

by 


T  (x)  •  min  T  . 

1  xeR  1 

A  level  a  test  of  HQ  verses  HA  is  given  by  reject  Hq  if 
T1  >  ‘l-o^l5* 

This  test  enjoys  all  the  same  properties  as  the  test  based 
on  T  described  in  Sections  2  and  3.  For  example,  the  test  has 
size  exactly  a  if  f^,  . ..,  f^  and  g  are  all  continuous  func¬ 
tions  and  there  is  a  value  of  8^  such  that  ffxJBj  ■  g(x)  for  one 
value  of  x  c  R  and  fj(x)8j  >  g(x)  for  all  other  x  e  R.  The  proofs 
of  these  properties  are  analogous  to  those  for  T  with  the  func¬ 
tion  fjCxjSj  -  f2(x)82  replaced  by  f^xJBj  -  g(x).  These  proofs 
are  not  given  herein. 


Table  1 


Power  of  the  Test  and  Percentage  of  Acceptances  and  Rejections 

1  2 
by  Shortcuts  for  Selected  Points  in  HQ. 


c 

I 

II 

III 

IV 

V 

0 

.0003 

.0003 

.0003 

.0003 

.0003 

99 

99 

99 

99 

99 

100 

100 

100 

100 

100 

.0027 

.0080 

.0037 

.0067 

.0087 

1 

100 

90 

98 

98 

99 

50 

33 

27 

50 

so 

.0053 

.0287 

.0110 

.0237 

.0327 

2 

100 

40 

97 

98 

99 

75 

22 

30 

39 

63 

.0060 

.0483 

.0267 

.0493 

.0523 

S 

100 

0 

97 

100 

100 

100 

16 

16 

32 

78 

.0060 

.0510 

.0460 

.0523 

.0523 

25 

100 

0 

99 

100 

100 

100 

15 

9 

30 

100 

.0060 

.0510 

.0523 

.0523 

.0523 

1000 

100 

0 

100 

100 

100 

100 

15 

S 

30 

100 

AFirst  (top)  entry:  estimated  power  ot  the  test 
Second  (middle)  entry:  percentage  of  acceptances  detected  by 

shortcut  in  Section  3.2.1 

Third  (bottom)  entry  :  percentage  of  rejections  detected  by 

shortcut  in  Section  3.2.2 


Column  I 

Column  II 

fj(x)B 

Column  III 

fjWe 

Column  IV 

fj(x)8 

Column  V 

fjWB 

f2(x)82  =  c(l  -  x2) 

*2  (*)®2  "  C*2 
f2(x)e2  -  c(x  ♦  l)2/4 

f2ooe2  a  c<x  *  1)^2 

f2(x)e2  »  C(-X2  ♦  2x  ♦  3)/4 
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Table  2a 

Power  of  the  Test  and  Percentage  of  Acceptances  and  Rejections 

1  2 
by  Shortcuts  for  Selected  Points  in 


c 

I 

II 

III 

IV 

V 

.0053 

.0290 

.0547 

.0180 

.0560 

5 

98 

94 

96 

99 

73 

38 

33 

51 

76 

21 

.0510 

.1460 

.2023 

.0990 

.2327 

1 

94 

87 

94 

99 

51 

45 

44 

63 

78 

36 

.5197 

.6897 

.7357 

.6030 

.7967 

2 

82 

81 

96 

99 

30 

65 

65 

81 

92 

59 

.9323 

.9647 

.9667 

.9387 

.9907 

3 

91 

93 

99 

100 

29 

90 

90 

97 

100 

88 

1.000 

1.000 

1.000 

1.000 

1.000 

S 

100 

100 

100 

100 

100 

*See  Table  1  footnote 


2 

Column  I 
Column  II 
Column  III 
Column  IV 
Column  V 


fjWBj  -  f2(x)S2  *  c 

fl(x)8l  '  f2(3t)S2  a  (x  *  1)2/4  *  c 

fl(x)Sl  '  f2(x)62  *  l_x2  *  2X  4  3)/4  *  C 
fjfxJSj  -  *2*x)e2  “  '*2  *  1  *  C 
fjtxiSj  -  f2(x)s2  *  x2  ♦  c 
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Table  2b 

Power  of  the  Test  and  Percentage  of  Acceptances  and  Rejections 

1  2 
by  Shortcuts  for  Selected  Points  in  H^. 


c 

I 

II 

III 

IV 

V 

.0053 

.0870 

.1453 

.0250 

.1510 

.5 

98 

94 

100 

100 

0 

38 

11 

79 

100 

17 

.0S10 

.2543 

.3167 

.1107 

.3680 

1 

94 

92 

100 

100 

0 

4S 

15 

87 

100 

28 

.  S197 

.7557 

.7700 

.6077 

.8460 

2 

82 

94 

100 

100 

0 

6S 

32 

96 

100 

57 

.9323 

.9677 

.9680 

.9387 

.9933 

3 

91 

99 

100 

100 

0 

90 

68 

100 

100 

88 

5 

1. 

0000  1.0000 

1.0000 

1.0000 

1.0000 

100 

100 

100 

100 

100 

^ee  Table 

1  footnote 

* 

^Column 

I 

-  f2«a2 

*  c 

5 

Column 

11 

fjWBj 

-  f2(x)e2 

*  (x  ♦  1)  «• 

2 

c 

Column 

III 

'  f2(x)02 

s  -x  ♦  2x  ♦ 

3  ♦  C 

Column 

IV 

-  *2^*2 

*  -4x‘‘  ♦  4  ♦ 

C 

Column 

V 

flWei 

- 

*  4xz  ♦  c 
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Figure  1 

Body  and  Kidney  Weight  for 
Healthy  (A)  and  Diabetic  (B)  Mice 
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Appendix 

Proofs  regarding  the  size  of  the  test  and  the  equivalence 
of  the  test  to  the  test  proposed  by  Tsutakawa  and  Hewett  (1978) 
and  Hewett  and  Lababidi  (1080)  are  given  in  this  appendix. 

The  size  of  the  test. 

Theorem  1:  Under  the  assumptions  of  our  model  the  test  has 
level  a,  i.c. , 


A.  1 


SUP  Pn  «  (T  >  ft(v))  5  a- 
(81*62)cH0  01,82  1 


Proof:  Fix  (6^,  e  HQ.  There  is  an  xQ  c  R  such  that 

£1^X0^81  5  f2^x0^62‘  Then*  with  Probability  one,  <  Q 

where 


fl(X0>bl  f2^X0^b2  '  <W81  ‘  f2^XQ^82^ 
S<^e(x0) 


Q  has  a  t  distribution  with  v  degrees  of  freedom.  So 


P81,e2(T  *  5  P«  "  ' 

Since  and  82  were  arbitrary,  A.l  is  true.  || 

Theorem  2:  Suppose  that  all  the  f^(x),  i  *  1,  2,  j  s  1,  •  ••»  P 

are  continuous  on  R.  If  there  exist  B^and  8^  sucb  *^at 

W81  "  W82  for  xo  c  R  80(1  >  f2(x)e2  £or 

all  other  x  e  R  then  the  test  has  size  exactly  a,  i.e.. 
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•UP  Pg  g 

(B1,fi2)€H0  V°2 


(T  > 


a. 


The  proof  of  Theorem  2  will  use  Lemma  1  which  can  be  proved 
using  standard  analysis  methods. 

Lemna  1:  Let  (x) ,  n  =  1 ,  2 ,  _ ,  be  continuous  functions  on 

a  compact  set  R.  Suppose  there  exists  an  xQ  e  R  such  that  gn(xQ) 
is  constant  (say  c)  for  all  n.  Suppose  gn(x)  increases  to  infinity 
as  n  +  «  for  all  x  *  xQ.  Then 


(A.2) 


lim  min  g  (x)  =  c. 
n-*“  xeR  n 


Proof  of  Theorem  2:  By  Theorem  1  it  suffices  to  show  there  exists 
a  sequence  (b"«  b") »  n  3  2,  ....  such  that,  (b",  b")  e  Hp 

for  n  ■  1 ,  2 ,  ....  and 


(A. 4)  lim  P  n  _n  (T  >  t  (v))  i  a  . 

n-w»  p1'b2  1-0 

The  estimates  b^,  i  «  1,  2,  can  be  written  as  b^  =  ♦  B^ 

where  Z^,  Z^  and  S  are  independent,  and  Z^  has  an  p .-variate 

normal  distribution  with  mean  0  and  variance-covariance  matrix 
2  -1 

a  D.  .  In  terms  of  these  quantities,  the  statistics  T  and  T 

X  a 

can  be  written  as 


T  -  T(Zj,  Z2,  S,  Bj,  B2)  =»  nin.Tx(Zl,  Z2,  S,  Bj, 

X€  R 


and 


,  f,(x)Z,  -  f,(x)Z  ♦  f  (x)B,  -  f-(x)B 

T  (Z  ,  Z  ,  S,  B.,  B,)  -  - ± - 1 - 2 - 1 - 1 - 2 - 2  — 

12 
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Consider  the  sequence  (Bj.  B2)  defined  by  g”  *  ng?  where  the 

g*  are  defined  in  the  statement  of  Theorem  2.  For  a  fixed 
1 

Pi  P- 

value  of  2j  £  R  ij  «  R"z  and  s  >  0,  define 
g^x)  =  Tx(zlt  z2,  s,  gj,  g2). 

The  g  (x)  satisfy  the  conditions  of  Lemma  1  since  1)  f . .  are 
n  13 

continuous,  2)  s/e(x)  >  0,  3)  fj(xQ)g"  =  f2^xo^2  an<* 

4)  fj(x)gj  -  f2(x)g”  increases  to  infinity  as  n  +  ®  for  all 

x  *  XQ.  By  Lemma  2,  lim  TfZj,  z^t  s,  g",  g”)  »  Tx  (Zj,  z2,  s,  g*,  g*) . 

Since  z^,  z 2»  and  s  were  arbitrary,  this  implies  that 

TfZj,  Z2,  S,  g",  g")  converges  to  Tx  (Zj,  Z2>  S,  g*,  g*)  with 

probability  one  and  hence  in  distribution.  Thus 

lim  P  n  „(T  >  t  (v))  »  lim  P(T(Z.,  Z  ,  S,  bJ,  Bp  >  ^ _o(v)) 
n-+°°  pl,p2  4  n-*®  4  4 

-  P(T  (Zx,  Z2.  S,  Bj,  Bj)  >  t10(v)) 

-a.  II 


Equivalence  with  tests  proposed  by  Tsutakawa  and  Hewett  (1978) 
sad  Hewett  end  Lababidi  (1980). 

j- 

Theorem  3:  Suppose  f, (x)B.  ■  g.„  ♦  I  6,  .x.  and  R  has  the 

— — — —  1  1  XU  lj  j 

from  R  *  (x:  x^*  s  x^  s  x?,  j  «  l,  ...,  k).  Consider  the  test 

which  rejects  Hrt  if  T*  >  t,  (v)  where  T*  ■  min  T  and  X*  is  the 

k  xeX* 

set  of  2  points  for  which  x^  is  either  x^*  or  x?.  Suppose 

a  s  .5.  Then  the  tests  based  on  T*  and  T  are  equivalent* 


Proof.  For  any  k  ♦  1  dimensional  vectors  and  b2  and  5  >  0, 
is  a  linear  function  of  (x^,  . ..,  x^)  divided  by  the  square 
root  of  a  quadratic  function  of  (x^,  x^)  which  is  posi- 

tive  for  all  (Xj,  ....  x^)  e  R  .  Such  a  function  has  the  pro¬ 
perty  that  T*  ■  min  T  £  0  implies  T*  =  min  T  a  T.  (This  is 
xeX*  x  xeR  X 

easily  proved  for  k  *  1  and  can  be  proven  for  general  k  by 

induction.).  For  any  bj,  b2  and  s  >  0,  T  s  T*  so  if  T  rejects 

Hq,  so  does  T*.  Suppose  b^,  t>2  and  s  are  such  that  T*  rejects 

H„.  Then  T*  >  t  (\j)  i  0,  since  a  S  .5,  so  T  =  T*  and  T 
0  l-o 

also  rejects  HQ.  |f 
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